Ferromagnetism in multi-band Kondo lattice model 
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The ferromagnetic spin exchange interaction between the itinerant electrons and localized mo- 
ments on a periodic lattice, studied within the so-called Kondo lattice model (KLM), is considered 
for multiband situation where the hopping integral is a matrix, in general. The modified RKKY 
theory, wherein one can map such a model onto an effective Heisenberg-like system, is extended to a 
multi-band case with finite bandwidth and hybridization on a simple cubic lattice. As an input for 
the evaluation of the effective exchange integrals, one requires the multi-band electronic self energy 
which is taken from an earlier proposed ansatz. Using the above procedure, we determine the mag- 
netic properties of the system like Curie temperature while calculating the chemical potential and 
magnetization within a self consistent scheme for various values of system parameters. The results 
are discussed in detail and the model is motivated in order to study the electronic, transport and 
magnetic properties of real materials like GdN. 

PACS numbers: 71.10.-w, 71.27. +a, 75.10.-b 



I. INTRODUCTION 

The original Kondo model^ with anti-ferromagnetic 
spin-exchange interaction between a single impurity spin 
in a non-magnetic background and the itinerant electrons 
of the host metal was used by Jun Kondo to explain the 
unusual temperature behaviour of the resistivity of the 
system. Its periodic extension with ferromagnetic ex- 
change interaction between a system of localized spins 
and a band of itinerant electrons is in the literature often 
referred to as the ferromagnetic Kondo lattice modefi^ 
or s-d££- or s-f model 6 . For the sake of uniformity, we 
ascribe it as the Kondo lattice model (KLM). Both the 
anti-ferromagnetic as well as the ferromagnetic alignment 
of itinerant and localized spins exhibits remarkable differ- 
ences in the physical properties of various real materials 
and has been a subject of intense theoretical studies in 
the past. 

For instance, the magnetic semiconductors (proto- 
types being the europium chalcogenides : EuX : 
X=0,S,Se,Te)£i£iii are known to have ferromagnetic ex- 
change coupling and demonstrate a spectacular temper- 
ature dependence of the band states. The redshift of the 
optical absorption edge in these materials upon cooling 
from T=T C to T=0 K is due to a corresponding shift 
of the lower conduction band edge&i -. A great deal of 
focus has been concentrated on studying the diluted mag- 
netic semiconductors with anti-ferromagnetic 1 ^ and fer- 
romangeti o 12 ' 13 exchange interaction with the purpose of 
achieving practical spintronic a 14 ' 15 applications. Apart 
from magnetic semiconductors, the local moment metals 
like Gd are known to have a ferromagnetic exchange 16 . 
But the exchange induced correlation and the tempera- 
ture dependent quasi-particle effects 1 -^ have lead to com- 
plex and hence controversial photoemission dat a 18 ' 19 . 
Other materials like the manganese oxides (mangan- 
ites) having pervoskite structures (the prototype being 
Ai_ x B x Mn0 3 where A=La,Pr,Nd and B=Sr,Ca,Ba,Pb) 
also have a strong ferromagnetic exchange interaction. 



They have a remarkable property called colossal magne- 
toresistance (CMR) 20 i 21 which enables them to dramati- 
cally change their electrical resistance in the presence of 
a magnetic field. Many theoretical models have been pro- 
posed in order to explain the existence of these effects. 
The earlier theoretical ideas were based up on the dou- 
ble exchange model 2 ^ which can be understood as one 
of the limiting case of Kondo lattice model (i.e. limit of 
strong Hunds coupling). Though recent theories^ 3 - have 
provided a step forward but its complete understanding 
is far from being explained by any current physical the- 
ories. As compared to the aforementioned compounds, 
the heavy fermion systems 2 ^ (mostly Ce compounds) are 
known to have anti-parallel alignment of the conduction 
electron and localized spins. They have been rigorously 
studied because of their extraordinary physical proper- 
ties 2 ^. 

In the above examples, the kinetic energy of the itiner- 
ant electrons is usually described within tight binding dis- 
persion of a single non-degenerate band i.e., single orbital 
atom per unit cell. But it is well known that the single 
band calculations are certainly not sufficient in order to 
have a complete understanding of unusual phenomenon 
in real materials 26 . One has to take into account the in- 
tra and inter- band interactions as well. 

The multi-band models are also of growing interest for 
exhibiting a wide range of phenomena like novel elec- 
tronic phases, magnetism and superconductivity 2 ^. For 
instance, it was found out using a two-band Hubbard 28 
model that there was a possibility of existence of ferro- 
magnetism around half-filling 29 in contrast to antiferro- 
magnetism in single band model 3 ^. The numerical stud- 
ies 3 ^ on ground state properties of multi-band periodic 
Anderson 3 ^ model revealed the minor role played by the 
competition between RKKY— and Kondo interactions 1 
again in contrast to the single band case^i. This moti- 
vates us to understand the physics behind the interplay 
between the kinetic and potential energy of multi-band 
Kondo lattice model and extend it for studying the elec- 
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tronic, transport and magnetic properties of real materi- 
als like GdN. 

This paper is organized as follows. In the following sec- 
tion we develop our theoretical multi-band model Hamil- 
tonian describing the physics due to the intra-atomic 
exchange interaction between the two sub-systems, i.e., 
itinerant electrons and localized spins, on a periodic lat- 
tice. In section Hi A[ we consider only the electronic part 
of the system while treating the magnetic part within 
molecular-field theory. Using an earlier proposed ansatz 
for multi-band self-energy^, we evaluate the electronic 
properties of interest like the density of states and band 
occupation number. In section lllBl we develop the mod- 
ified RKKY theory^ wherein we integrate out the charge 
degrees of freedom of the itinerant electrons thereby 
mapping the multi-band Kondo lattice Hamiltonian onto 
an effective Heisenberg- like spin Hamiltonian. In sec- 
tion |TTT1 using the procedure described in earlier section 
we determine the magnetic properties of the system like 
Curie temperature while calculating the chemical poten- 
tial and magnetization within a self consistent scheme. 
We discuss the results as obtained for various values of 
system parameters. In section IIV1 we summarize and 
conclude our findings. 

II. MODEL AND METHODS 

A. Electronic sub-system 

In this section, we present a brief description of the 
theoretical model used in our calculations. The details 
of the many body analysis along with a model calculation 
and limiting cases are explained elsewhere^. The multi- 
band Kondo lattice model (KLM) Hamiltonian mainly 
consists of two parts ; 

H = Hkin + Hi n t (1) 

where 

H Hn = T ilAa.o c i&° ( 2 ) 

and 

Hint = ~"^"E/^* ' <T ) cr ' ' C Lcr' C iao- (3) 

ia 
tr'a 

Hkin denotes the kinetic energy of the itinerant elec- 
trons with T^ 3 being the hopping term which is con- 
nected by Fourier transformation to the free Bloch ener- 
gies e a/3 (k) 

T n = ^E^( k ) e-^-^) (4) 

k 

while c\ aa and Ci aa are the fermionic creation and annihi- 
lation operators, respectively, at lattice site Ri. The latin 



letters (i,j,...) symbolize the crystal lattice indices while 
the band indices are depicted in Greek letters (a,f3,..) 
and the spin is denoted as a(= f, I)- 

H^t is an intra-atomic exchange interaction term i.e., 
a local interaction between electron spin <x and local mo- 
ment spin Si. Using second quantization for electron spin 
{n-iaa — cl aa Ci aa ), the interaction term is further being 
further split into two subterms. 

Hint — ~ 7T y ] ( z <rSj c iaa C ia& + &i c ia-a C iacr) (5) 

The first describes the Ising type interaction between 
the z-component of the localized and itinerant carrier 
spins while the other comprises spin exchange processes 
which are responsible for many of the KLM properties. 
J is the exchange coupling strength which we assume to 
be k-independent and Sf refers to the localized spin at 
site Ri 

S° = S*+tz a Sf ; z T =+M i = -l (6) 

The Hamiltonian in eq. |T]) provokes a nontrivial many 
body problem that cannot be solved exactly. Using the 
equation of motion method for the double-time retarded 
Green functional 

GZAE) = {{ci^,cl va )) E (7) 

where l,m and [i^v are the lattice and band indices re- 
spectively, we obtain higher order Green functions which 
prevent the direct solution. Approximations must be con- 
sidered. But a rather formal solution can be stated as 

G kCT (£) - \{E + iO+)T -?(k) - tuAET 1 (8) 

where for simplicity we exclude the band indices by rep- 
resenting the terms in a generalized matrix form on sym- 
bolizing a hat over it 

G lma (E) = jjY,GUE) e-^< R '- R ^ (9) 

k 

The terms in eq. (|8]) are explained as follows : / is an 
identity matrix. + is small imaginary part and e(k) is 
a hopping matrix with the diagonal terms of the matrix 
exemplifying the intra-band hopping and the off-diagonal 
terms denoting the inter-band hopping. The self energy, 
Sko-IE), containing all the influences of the different in- 
teractions being of fundamental importance, can be un- 
derstood using site representation : 

(([H inU cv,]_; cL CT » = E ^{E)G^ na {E) (10) 

pi 

Now we are left with a problem of finding a multi-band 
self energy ansatz in order to compute the Green function 
matrix and thereby calculate the physical quantities of 
interest like the quasi-particle spectral density (SD) 

S ka (E) = --ImTr{G^{E)) (11) 
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and the quasi-particle density of states (Q-DOS) 



Pa(E) 



Nh 



which would yield the band occupation number 

/oo 
dEf-{E) Pa {E) 
-OO 



(12) 



(13) 



where 



f-(E) 



1 



is the Fermi function and /i being the chemical potential 
or the Fermi edge. 

According to our many body theoretical analysis^, we 
utilize the multi-band interpolating self energy ansatz 
(ISA) which is well defined in the low carrier density 
regime 3 ^, for all coupling strengths and satisfying one 
limiting case of the model namely that of ferromagnet- 
ically saturated semiconductor. The ansatz is given as: 



E CT (£) = J -M_ a I+ ^a_ a G- a (e + 3 -M_ 



r-io. 



E + J -M_ c 



(14a) 



where 



M a = z a {S z ); 



S{S+l)+M a (M a + l). (14b) 



and the bare Green function matrix is defined as 
G CT (£) = ^£[(£ + iO + )7-e(k)]- 1 



The first term in Eq. (|14ap , which is exact in the weak 
coupling limit, represents an induced Stoner splitting of 
the energy band proportional to the / spin magnetization 
(S z ). The second term is dominated by the consequences 
of spin exchange processes between itinerant electrons 
and localized / moments. 

As seen in Eq. (|14a[) . we are interested only in the local 
self energy 



determine the influence of inter- band exchange on the 
conduction band states so as to study the electronic cor- 
relations effect o 35 ' 38 . But in order to study the effect of 
itinerant electron subsystem on the localized subsystem 
and vice versa, we need to calculate the magnetization 
within a self-consistent manner as shown in the following 
section. 



B. Magnetic sub-system 



In Section QTA] we did not consider a direct exchange 
interaction between the localized / spins. But if one is in- 
terested in determining the magnetic properties of multi- 
band Kondo lattice model then both the sub-systems (lo- 
calized as well as itinerant) should be solved within a self 
consistent scheme. Therefore, we would like to take into 
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Figure 1: (Color online) An effective indirect exchange, J^ ff , 
between localized / spins (red arrows) mediated by intra- 
atomic exchange, J, due to itinerant electrons. Ef denotes 
the Fermi edge. 

account an effective indirect coupling, J?j ff , between the 
localized / spins and the itinerant electrons within the 
so-called modified RKKY 36 formalism as shown in Fig- 
ure [1] 

Consider the multi-band Kondo lattice Hamiltonian, 
i.e. Eq. ([TJ) , which can be written in the following equiv- 
alent form, 

H = Hkin + Hi n t 



j 

2N 



ia kq 

cr'er 



while the wave-vector dependence of the self energy is 
mainly due to the magnon energies appearing at finite 
temperature. In order to evaluate only the itinerant 
electron subsystem, we can neglect this wave-vector de- 
pendence. The localized magnetization (S z ) can then 
be considered as an external parameter being responsi- 
ble for the induced temperature dependence of the band 
states. Thus in a non-self consistent way, it is possible to 



where all the terminologies remain the same as explained 
in Section III Al The components of the band electron 
spin operator <r are the Pauli spin matrices. 

The main idea of the modified RKKY theory is to 
transform the above Kondo-like exchange Hamiltonian of 
the conduction electrons into an effective Heisenberg-like 
spin exchange Hamiltonian of the / spins by averaging 
Hi n t in the subspace of the conduction electrons 
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(( )): 

{Hint) = Hf 



J 

"2N 



ia kq 
a' a 



(15) 



This is achieved by closely following the treatment as 
given in Ref . l36l. For averaging procedure the / spin oper- 
ators are to be considered as c numbers. The expectation 
values ( ) in Eq. (|15|) may still have operator proper- 
ties in the / spin subspace and therefore do not vanish 
for q ^ and a ^ a' . We would like to obtain ( } via 
the spectral theorem with the help of appropriate Green 
function as given below : 



G k!k+q(^) = (( C k QCT ;c3 c +q/3 CT '))£ 



(16) 



Its equation of motion can be obtained in the usual 
way^i and is given as 



(17) 



The above equation can be iterated up to any desired 
accuracy producing spin products of the type : 
(Sj • <J)oo" , (Sj • <t) ct " ct '" , (Si ■ cr) a "i a "" 
On excluding the band indices in Eq. (|17[) by represent- 
ing the terms in a generalized matrix form on symbolizing 
a hat over it we get 

EG^(E) = <5 k , k+q <W/ + ^k)G^ k+q (£) 



2N 



ik'cx'' 



(18) 



Rearranging the terms in Eq. (|18p yields 
[ET— e(k)] G k ^ +q (E) = <5 k ,k+ q <W/ 

E e-^-^(S i .crU,G^ +(i (E) 

(19) 

For symmetry reasons, we write down the equation of 
motion for G k k °^ q (i?) in an alternative way where the 

second operator , c k+q/3cr , , in Eq. fTo]) is the "active" 
operator.: 

EG& +(l (E) = 5 k , k+q <W?+G^ k+q (i?)e(k+q) 

~W E e-^'-^^.a^G^E) 

(20) 



ik' ' a n 



and again upon rearranging the terms in the above equa- 
tion we get 

G k 7 k+q (£)[£7-e(k+q)] = <5 k>k+q <W7 



-E 

O.N 



-i k'-(k+q) -Ri 



2N 



(Si • cr) (T // (T /G' k ff k ,(7;) 



ik'cr" 



(21) 



Now let us define the following Green function of the 
'free" electron system 



[eT-?m] = (g£\e))- s 



(22) 



(23) 



J 

"2N 



[EI-eik+^^iG^E))- 1 
Subsituting Eq. (J22J) in Eq. (19]) gives 
(g£\e)) 1 G k CT k+q (£') = <5 k , k +q<W7 

(24) 

while substituting Eq. (J23J) in Eq. (J5TJ) yields 

G k CT k+q (S)(G| t j. q (-B)) 1 = <5 k , k+q <W? 

E e -< k '-( k +i))- R '(S, • o-V^G^i?) 

(25) 

Now upon multiplying G k O) (£0 from left to Eq. (24|) 
we get 

G k CT k +q(E) = #k,k+q<WG k ° \E) 

1 E ^ l(k - k ' ) - R '(s l ^) CTCT "G| c 0) (i : ;)G^ +q ( J B) 



J 

"2iV 



27V 



(26) 



And on multiplying G k °|_ (£) from right to Eq. ([25]) 
we obtain 

G k CT k+q (£) = <$k,k+q<WG k °[ q (.E) 

E e-< k '- (k+q) )^(S i . CT ) CT ^,G^ k :(75)G k °j q ( J B) 

(27) 



Let us make the following crucial first order approxi- 
mations for the Green functions 



Gu%'(E) ~ 5 CTCT "5 k , k 'G k(T (7^) 



(28) 



where 



Gk'M+q( E ) ~ ^-"cr'fSk'.k+qGk+qo-'GE) (29) 



G kCT (£) = [£/ - e(k) - £ CT (£!)] 1 (30) 
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G k+qff ,(£) = [EI- e(k+q) - t a ,(E)\ 



(31) 



The renormalization by the interacting Green functions 
as performed in Eq. (f2"B")) and Eq. (j2"9")) should be a sensible 
approximation since it is observed that if those interact- 
ing Green functions are replaced by the free Green func- 
tions, Eq. (|22|) and Eq. ((23|) respectively, then it leads to 
the correct low- J (i.e. RKKY) behaviour. On substitut- 
ing Eq. (29I in Eq. ([26]) we obtain 



Hf = 

J 

2nN 



/OO 
dEf-(E)Y,G { £\E) 



J 



8irN 2 



E 

ijqcrcr' 



G 



k,k+q 



ImTr r dEf_(E)^A^ +(l (E) 

J -oo k 

(37) 



J 



■ £ e^ R * (Si ■ <r) aa , (£)G k+qCT , (£) 

(32) 

while substituting Eq. (J28J) in Eq. (|27|) gives 



i.e. 



Gkjk+q^) - £q,o<WG k ° (E) 

J ^ e ^(S. t ■ <r)^G kCT (i?)G^ q (£;) 



-ijqcrcr' 



q 

(38) 



2N 



On adding Eq. 



(33) 



and Eq. (f3"3"|) we get 



G^k+qt^) - ^q.O^'G^-E) 

J .^ e ^(S. 1 . ( T) ffCT ^ k ff k+q ( J B) 



AN 



where 
(4 0) 
and 



1 f°° 1 

) = —IrnTr (£) - £ G k 0) (£) (39) 

^ — OO i. 



1 f 00 1 / 

■-/mTr/ d£/_(£)-^A k CT k+q (£) (40) 



(34) 



where 



If we perform spin summations on r.h.s of Eq. (|38[) we 
obtain 

J 



^k a k +q (£) = (G { °\E)G W (E) + G kCT (£)G k °[ q (£;)) Hf = «n< 0) ) - <n< 0) »S ; 



(35) 

For the effective spin Hamiltonian in Eq. (fTS"|) we need 
the expectation value (c k _|_ qQ , ./Ck QC r), which we express in 
terms of the trace of imaginary part of the Green function 
Eq. (f34|) by exploiting the spectral theorem^! : 



8iV 



V^g-iq^Ri-Rj) 



jy E^ C k+qc«7' C kaCT 
k 

1 f°° 

-—ImTr / dEf- (E) £ G£ ff k+q (i?) 
= <5 q ,o<W (^j ImTr J dEf- (E) £ G k ° } (E) 



D^S~S+ + D^S+SJ + (X)TT + DU)5?5; 



(41) 



where the spin operators, Sf (£ = +,-,z ), satisfy the 
usual commutation relations. The first term in the above 
equation is exactly zero since the free system is un- 
polarizcd. This finally yields an effective anisotropic 
Heisenberg- like spin Hamiltonian which can be written 
as follows 



J 

' 4ttN 2 



E 



e iq - Ri (Si-<7)<r<r' * 

/CO 
dEf^E^Ki+^E) 
-oo k 

(36) 

On substituing Eq. (f3"6")) in Eq. (fT5|) we get 



r(l) c— c+ I 7l z i q+ c— I tW qz qz 



1+ _L 7-( 2 ) C+ < 



7(3), 



(42) 



where 



^-^E^^ ' " V{H " 1 (» = 1A3) 

(43) 



with 



Using Eqs. fl49|), (|50|). (|5T|) and (J52J) we obtain 



^ (1) (q) = -y^ 
J< 2 >(q) = -£^ T 

j(3) (q ) = -^( I) n + jD u ) 



(44) 
(45) 
(46) 



are the effective exchange integrals which via G^ a are 
functionals of the conduction electron self-energy thereby 
getting a temperature and carrier concentration depen- 
dence. In order to obtain effective isotropic Heisenberg- 
like spin Hamiltonian, one can prov e22- that = 

and + = 2D^ which will result in jW (q) = 
J( 2 )(q) = j( 3 )(q)/2 = J off (q). We finally get 



- ( Si Sf + SfSj 



QZ QZ 



(47) 



Now so as to deterimine the / spin magnetization we 
follow along the lines of Calient that results in 



_ (S- y )(l + y ) M " + (S + l + y )^ 



(1 + *) 



25+1 _ (p2S+l 



where if(S) can be interpreted as average magnon num- 
ber 



i 



N ^ e E(q)/k B T _ 1 

q 



(49) 



depending on S via magnon energies E(q) which can be 
obtained using the spin Green function^ and is given by 



E(q) = 2{S%J cS (0)~ J cff (q)] 



(50) 



with J off (0) = J cff (q = 0). In section |HA1 we observed 
that the magnetization ((5 Z )) appears in electronic self 
energy. And this self energy is used to calculate the ex- 
change integrals which along with {S z } enter the magnon 
energies. These magnon energies in turn also appear in 
(S z ). Thus we have found a closed system of equations 
that can be solved self consistently for all quantities of 
interest, in particular those which tell us about the mu- 
tual influence of electronic and magnetic properties of 
the exchange coupled system of itinerant electrons and 
localized / spins. 

One of the central quantities in magnetic system is the 
Curie temperature, T c , which can be ascribed to the tem- 
perature for which (S z ) — ► 0. On expanding Eq. (|4"5|) in 



¥>(S) 



we get 



(S z ) = 



5(5 + 1) 



For (S z 



MS) 

we have 

e E(q)/k B T _ 1 



o 



1 



MS)) 2 
E(d) 



k B T 



(51) 



(52) 



25(5 + 1) 
3k B 



N ^ 



1 



[jeff(0)-Jeff(q)]y Ti 



(53) 

We can evaluate Eq. (f53|) within a self consistent 




GAE), £^(-E), GW-E), pAE) 




Figure 2: Flowchart exhibiting the self consistent determi- 
nation of Curie temperature, T c . The terminologies are as 
explained in the text. 

cycle as shown in Figure [2] which can be understood 
as follows. In our analysis we consider a two band 
model (a and v — 1,2) which can be generalized to a 
n-band model. The single particle energies e(k) is then 
represented in a 2 x 2 matrix where the diagonal and 
non-diagonal terms are considered to have the following 
form. e n (k) = — ^-(cos(k x a) + cos(k y a) + cos(k z a)) 
, e 12 (k) = e 21 (k) = V (local hybridization, LH) or 
e 12 (k) = e 21 (k) = Ve 11 (k) (non-local hybridization, 
NLH) and e 22 (k) = E + e n ( k )- A1 «ng with the 
single particle energies, the intra-atomic exchange J, 
quantum spin number S and bandwidth W act as input 
parameters in order to evaluate the free propagator, 
self-energy and full propagator. Then for a particular 
band occupation, n, and an initial temperature T™, 
the Fermi edge, y, which yields the correct value of n, 
is determined self consistently. Thus, fixing upon the 
Fermi edge one evaluates the temperature dependent 
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exchange integrals which gives the T c through Eq. ([53)1 . 
If the obtained temperature is within convergence limit, 
e, then it is the resulting T c for particular J and n. 



III. RESULTS & DISCUSSION 

On optimizing the numerical factor—, we evaluate the 
Curie temperature, i.e. Eq. (I53|) . for various configura- 
tions of model parameters (J, n and V). We consider 
the spin quantum number S = | (localized moment of 
Eu 2+ , Gd 3+ ) and the center of gravity of the second 
band shifted by amount, Eq — 0.25 eV. First we try to 
reproduce the single band result for different values of 
band occupation as was previously obtained but using 
another electronic self-energy^. This will give us some 
confidence on the working of the algorithm. 

Figure [3] shows the dependence of T c on the strength 



local hybridization 



non-local hybridization 





Figure 3: (Color online) The dependence of Curie tempera- 
ture (T c ) on intra-atomic exchange (J) for different values of 
band occupation. The exhibited results are for single and two 
band KLM with local (left panel) and non-local hybridization 
(right panel) on a simple cubic lattice. 

of intra-atomic exchange, J, for single and two band 
KLM. The left panel describes the results with local 
hybridization (LH) while the right panel for non-local 
hybridization (NLH) on a simple cubic lattice. The 
bandwidth of both the bands, W, is taken to be 1.0 eV. 
The calculations are carried out for different values of 
band occupation, n, and hybridization, V. 

Though both the graphs look quite similar but there 
are marked differences especially in the limit of strong 
coupling and low band occupation. We first discuss the 
general behavior. It is observed that initially the T c rises 
sharply with increasing J. For weak coupling and small 



band occupation (n=0.05), the usual RKKY mechanism 
is observed. However for higher band occupation, J 
is observed to exceed a critical value in order to allow 
ferromagnetism. Furthermore, with increasing J the 
critical temperature is observed to be deviating more and 
more from the RKKY behavior (i.e., long- range order) 
and finally reaches a saturation. The calculations done 
within single band model are comparable with previous 
calculations^- obtained using a different self-energy. 

As shown in Figure [51 the results for T c in case of 
two unhybridized bands (V=0.0) are in comparison with 
that of the one band situation for low band occupation 
due to similar low energy paramagnetic density of states 
at the Fermi edge. And with increasing band occupation 
the results in both the situation (local and non-local 
hybridization) differ drastically. The values of Curie 
temperatures are found to be higher and increasing with 
increasing hybridization strength and band occupation 
for two hybridized bands as compared to the unhy- 
bridized or one band model. It can be understood as 
follows. 

Figure 2] shows the low energy window of the 



J,=0.20eV,<S>=0.0, LH 




NLH 


1 

n=0.05 


1 1 1 1 1 1 
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Figure 4: (Color online) Lower edge of paramagnetic density 
of states of two band KLM with LH (upper panel) and NLH 
(lower panel) on a sc lattice for J— 0.20 eV. The curves and 
vertical lines in yellow, orange and red are for V=0.0, 0.1 and 
0.2 respectively. The vertical lines denotes the Fermi edge 
with dashed-dotted and dashed lines for a band occupation, 
n, of 0.05 and 0.30 respectively. 

paramagnetic density of states of two band KLM with 
local and non-local hybridization shown in upper and 
lower panel respectively. The calculations are performed 
for J=0.20 eV. The curves and vertical lines in yellow, 
orange and red are for V=0.0, 0.1 and 0.2. The vertical 
lines signifies the Fermi edge with dashed-dotted and 
dashed lines for a band occupation of 0.05 and 0.30 
respectively. The lines within the circle represent the 
Fermi edge for different values of hybridization but for 
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the same value of band occupation. 

It is to be noted that in general the shape of the 
density of states is dependent on the band occupation 
which is a consequence of electronic correlation effects. 
But the self-energy which we consider in our calculations 
is independent of band occupation or rather well defined 
only in the limit of low band occupation. Since we 
restrict ourselves to this limit, thus we have density of 
states dependent only on the strength of hybridization. 
The band filling is determined by the placement of 
Fermi edge which is obtained self-consistently in our 
calculations. In order to improve over the restricted 
limit, we can consider the band occupation dependent 
self-energy as given in Ref. But this is not the aim 
of the present paper. 

We observe that for low band occupation, the density 
of states at the Fermi edge are slightly different for 
different values of hybridization. But with increasing 
band occupation, the density of states at the Fermi 
edge changes abruptly for local as well as non-local 
hybridization. 

Though for high band occupation the values of T c 



n=0.05 n=0.10 n=0.30 n=0.50 




0.4 0.8 0.4 0.8 0.4 0.8 0.4 0.8 



J (eV) -> 

Figure 5: (Color online) The same as in Figure[3]but only for 
LH and an additional result for band occupation of n=0.50. 

in case of two hybridized bands are higher as compared 
to unhybridized or one band situation but for low band 
occupation and strong coupling limit, an interesting 
feature is observed as shown (encircled) in Figure [5] 
Since it is only observed in case of two locally hybridized 
bands so we do not consider the case of non-local 
hybridization. It is noted that in the limit of strong 
coupling the T c starts decreasing for two hybridized 
band system as compared to two unhybridized band 
model. Even though increasing hybridization increases 
the bandwidth and therefore the kinetic energy of the 
itinerant electrons but it is only effective for higher band 
occupation. In the regime of low band occupation and 
strong coupling the short range order due to strong 



intra-atomic exchange and local hybridization decreases 
the kinetic energy leading to localization of the carrier. 
This results in decrease of the paramagnetic density 
of states at the Fermi edge thereby reducing the T c 
to the value of the one band case which is encircled 
in the left-most panel in Figure [5] But as mentioned 
earlier, upon increasing the band occupation (moving to 
the right panels) we observe that the T c increases with 
increase in hybridization (bandwidth^) and also in the 
limit of strong coupling due to the presence of more 
delocalized electrons. It can be again understood within 
the picture of the density of states. 

Figure [6] shows the low energy spectrum of para- 

J st ^».80eV, <S*>=0.0, LH 
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\( ( , 




1 
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Figure 6: (Color online) Lower edge of paramagnetic density 
of states for two locally hybridized band KLM on a sc lattice 
for J=0.80 eV. The curves and vertical lines in yellow, orange 
and red are for V=0.0, 0.1 and 0.2 respectively. The vertical 
lines denote the Fermi edge with dashed-dotted and dashed 
lines for a band occupation, n, of 0.10 and 0.50 respectively. 

magnetic density of states for two locally hybridized 
band KLM on a simple cubic lattice for J=0.80 eV. The 
curves and vertical lines in yellow, orange and red are for 
V=0.0, 0.1 and 0.2. The vertical lines denote the Fermi 
edge with dashed-dotted and dashed lines for a band 
occupation of 0.10 and 0.50 respectively. As observed 
for n=0.10 that the density of states at the Fermi edge 
are slightly different giving rise to marginal difference in 
T c with increase in hybridization. But with increasing 
band occupation (n = 0.50), the density of states at the 
Fermi edge differ drastically giving rise to quite different 
Curie temperatures. 

Another feature which we observe in the strong 
coupling limit is that the density of states tend to 
seperate out with increasing value of hybridization as 
shown in Figure [6l But the change from V=0.1 to 
V=0.2 is quite sudden. In order to have a close look at 
it, we plot the paramagnetic density of states for J=0.80 
eV and intermediate values of hybridization from V=0.1 
to V—0.2 as shown in Figure [7] It is noted that there is 
an increase in bandwidth with increasing hybridization. 

Further interesting characteristic is observed for band 
occupation of n=0.50. It is seen that the trend of higher 
and increasing value of T c for two locally hybridized 
band model is reversed as shown in the right most panel 
of Figure As the Fermi edge keeps on moving to 



9 



— v^o.oooo 

V^O.IOOO 

V=0.1500 

— V-0.1750 
V-0.1S75 

v-y.2«oo 




values of hybridization. We observe the long- range 
RKKY kind of oscillations 33 for weak coupling in case 
of local and non-local hybridization. In case of strong 
coupling, the local short range order is more strong. 
In that case, the exchange integrals get converged very 
quickly within a short distance. 



IV. SUMMARY & CONCLUSION 



Figure 7: (Color online) Lower edge of paramagnetic density 
of states for two locally hybridized band KLM on a sc lattice 
for J=0.80 eV and for different values of hybridization. 



higher energies with an increase in band occupation, 
the paramagnetic density of states at the Fermi edge 
keeps on changing. This results in an observed change 
in the T c . A similar pattern of decrease in the value 
of Curie temperature with increase in band occupation 
and strength of hybridization is also observed in case 
of non-local hybridization. The explanation lies similar 
to what we have explained earlier in case of local 
hybridization. 

One the other hand since T c is directly related to 
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Figure 8: (Color online) The indirect exchange integrals 
shown as a function of distance for local (upper panel) and 
non-local hybridization (lower panel). The shown results are 
for two different values of intra-atomic exchange and band 
occupation but for three different values of hybridization. 

effective exchange integrals, J° ff , it is also interesting 
to notice the behavior of these exchange integrals for 
different values of n, J and V and for local as well as 
non-local hybridization. Figure [5] shows the dependence 
of indirect effective exchange integrals on two different 
parameter configuration of J, n and for three different 



In this paper, we studied the magnetic properties of 
the multi-band Kondo lattice model Hamiltonian which 
describes the intra-atomic exchange interaction between 
itinerant electrons and localized spins on a periodic 
lattice. In section III A[ we considered only the elec- 
tronic part of the system. Using an earlier proposed 
ansatz for multi-band self-energy 3 -^, we evaluated the 
electronic properties of interest like the density of states 
and band occupation number. In section IIIB1 we de- 
veloped the modified RKKY theory 3 ^ wherein we inte- 
grated out the charge degrees of freedom of the itinerant 
electrons thereby mapping the multi-band Kondo lattice 
model Hamiltonian onto an effective Heisenberg- like spin 
Hamiltonian. In section HTT1 using this procedure we de- 
termined the magnetic properties of the system like Curie 
temperature (within Random Phase Approximation) for 
various values of system parameters while calculating the 
chemical potential and magnetization within a self con- 
sistent scheme. 

We found that the T c as a function of intra-atomic ex- 
change J for a two band KLM remains qualitatively the 
same for local as well as non-local hybridization between 
both the bands except for the limit of low band occu- 
pation and strong coupling. For higher band occupation 
and increase in strenght of coupling as well as hybridiza- 
tion we find that T c increases until the band occupation 
of n=0.5 from where the trend is reversed. All these can 
be explained using the paramagnetic density of states and 
its behaviour at the Fermi edge. It is mainly due to the 
interplay between kinetic and potential energy. In case 
of strong coupling, the T c is oscillating in its dependence 
on the band occupation. 

Such an analysis can be very useful in order to under- 
stand the physical properties of real materials described 
within the multi-band models like the manganites. These 
materials are known to have a strong intra- atomic ex- 
change coupling behavior. Or the analysis can be equally 
handful for the rare- earth metals which are known to 
be described within the weak or inter- mediate intra- 
atomic coupling regime. It would be equally encouraging 
to carry out the similar investigation for two different 
bandwidths of both the bands since the correlation ef- 
fects scale as where W being the bandwidth. 

We would also like to apply 4 ^ the multi- band modified 
RKKY theory in order to understand the basic mecha- 
nism behind the observed ferromagnetism in Gd N 45 i 46 . 
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